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This work presents the results of an experimental study on supercooled droplet 
deformation and breakup near the leading edge of an airfoil. The results are compared to 
prior room temperature droplet deformation results to explore the effects of droplet 
supercooling. The experiments were conductedin the Adverse Environment Rotor Test Stand 
(AERTS ) at The Pennsylvania S tate University. An airfoil model placed at the end of the rotor 
blades mounted onto the hub in the AERTS chamber was moved at speeds ranging between 
50 and 80 m/sec. The temperature of the chamber was set at -20°C. A monotonic droplet 
generator was usedto produce droplets that fell from above, perpendicular to the path of the 
airfoil. The supercooled state of the droplets was determined by measurement of the 
temperature of the drops at various locations below the droplet generator exit. A temperature 
prediction code was also used to estimate the temperature of the droplets based on vertical 
velocity and the distance traveled by droplets from the droplet generator to the airfoil 
stagnation line. High speed imaging was employed to observe the interaction between the 
droplets and the airfoil. The high speed imaging provided droplet deformation information as 
the droplet approached the airfoil near the stagnation line. A tracking software program was 
used to measure the horizontal and vertical displacement of the droplet against time. It was 
demons trated that to compare the effects of water supercooling on droplet deformation, the 
ratio of the slip velocity and the initial droplet velocity must be equal. A case with equal slip 
velocity to initial velocity ratios was selected for room temperature and supercooled droplet 
conditions. The airfoil velocity was 60 m/s andthe slip velocity for both sets of data was 40 m/s. 
In these cases, the deformation of the weakly supercooledand warm droplets did not present 
different trends. The similar behavior for both environmental conditions indicates that water 
supercooling has no effect on particle deformation for the limited range of the weak 
supercooling of the droplets tested and the selected impact velocity. The assumption of a 
cons tant surface tension value was further supported by the equal trend of the Bond number 
obtained for supercooledand room temperature droplets. 
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Nomenclature 


avi = Audio video interleave movie format 

DBKUP 002 = Designation for Airfoilused in the experiment 

INTA = Instituto Nacionalde Técnica Aeroespacial 

LED = Light Emitting Diode 

Slip Velocity = Relative velocity between thedroplet and the air: (Varoptet - Vair) 


I. Introduction 


Aircraft manufacturers have reported in-flight observation of droplet breakup near the wing surfaces of large 
transport aircraft’. This observation is important for in-flight icing because if large droplets break up before impinging 
on aircraft surfaces, the smaller droplets thatresult may not impinge on the airfoil potentially affecting the ice accretion 
process. If droplet breakupis a factor in the ice accretion process, currentice accretion codes willneed to be modified 
to account for the breakup phenomena. 


In 2005 the Federal Aviation Administration (FAA ) sponsored a computational study by Wichita State University 
(Tan, Papadakis and Sampath’) to assess the effects of aerodynamic forces on water droplets near the leading edge of 
an airfoil and the leading edges of the slat and flap elements of a high-lift airfoil. Two-dimensional numerical 
computations of droplet breakup were done using the TAB model in the commercial code FLUENT. The TAB model 
used included only vibrational and bag models of breakup. The effects of chord size and droplet diameter were s tudied 
fora NACA 0012 geometry forchords of 0.91 m (3 feet) and 6.1 m (20 feet) and droplet sizes of 100, 500 and 1000 
micrometers. For the case of a 6.1-m (20-foot) chord airfoil, droplet breakup was studied for a three element airfoll 
high lift systemin a landing configuration. The study indicated that droplets may breakup in regions with severe 
pressure gradient. For the 0.91-m (3-foot) chord single element airfoil the droplet breakup occurred aft of the leading 
edge. For the 6.1-m (20-foot) chord single element airfoil the droplet breakup occurred near the stagnation region 
only. Large droplets were found to be more susceptible to breakup than smaller ones. They found that when droplet 
breakup occurred near the airfoil surface there was insufficient dis tance between the airfoil wall and the location where 
the droplet breakup was initiated for the droplets to achieve complete breakup before hitting the airfoil. The trajectory 
followed by droplets larger than 500 micrometers were “ballistic” in nature. The study pointed out the lack of 
experimental data in this area and recommended experimental tests to assess the effects on droplet breakup and 
breakup modes fromthe pressure gradient and therelative droplet-gas velocity near the leading edge of an airfoil. 


In 2007 the National Aeronautics and Space Administration (NASA) Glenn Research Center and the Instituto 
Nacional de Técnica Aeroespacial (INTA) in Madrid, Spain, began an experimental research programto obtain dropket 
breakup data on an airfoil configuration. A droplet breakup rotating rig? was designed and built at the INTA 
installations near Madrid. The first sets of experiments were conducted at low speeds (15-66 m/s ) in the fall of 2008. 
Using the experience gained fromthe low speed experiment thetestrig was modified to attain speeds upto 90 m/sec. 
Additional experiments were conducted between 2008 and 2012 for velocities from 50 to 90 meters per second and 
droplet sizes from 100 to 1800 micrometers’. The results allowed the gain of understanding on the phenomena and 
the behavior of the important parameters. All the experiments were conducted at roomtemperature. Since in in-flight 
icing the droplets are supercooled, the effect of supercooling was not answered at the time and remained an open 
research question. 


The present work reports the results of a prelimmary experiment conducted in the Adverse Environment Rotor 
Test Stand (AERTS) facility at The Pennsylvania State University to determine the effect of supercooling on dropkt 
deformation and breakup as they approach the leading edge of an airfoil. In the AERTS facility, air temperatures 
below freezing can be reached and the droplets can be supercooled. An airfoil model placed at the end of the rotating 
arm was moved at speeds of 50 to 80 m/sec. A monosize droplet generator produced droplets from 200 to 700 
micrometers that were allowed to fall from above, perpendicular to the path of the model airfoil at a given location. 
The temperature of the facility was maintained at -20°C. The droplets were confirmed to be supercooled based on the 
measurement of the temperature of droplets at various distances fromthe exit of the droplet generator using an infrared 
camera. A droplet temperature prediction code was alsoused to verify the supercooling of the droplets based on the 
environmental conditions, water temperature at the droplet generator, and exit velocity. The prediction tool also 
confirmed the droplets could be supercooled at the impact location given the conditions provided. High s peed imaging 
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was employed to observe the droplet deformation and breakup. A tracking software program was used to measure, 
from the high speed movies, the droplet horizontal and vertical dis placement against time. The horizontal dis placement 
data was curve fitted to obtain the velocity and acceleration. The velocity and acceleration, together with experimental 
values of the air velocity at the locations of the droplet, were used to calculate the W eber and Reynolds numbers along 
the path of the droplet from beginning of deformation to breakup and/or hitting the airfoil. The parameter behavior 
and the droplet deformation were compared to previous results fromstudies done at roomtemperature. 


In the present work, for the first time, droplet deformation and breakup has been conducted using supercooled 
droplets approaching the leading edge of an airfoil. Although prelimmary in nature given the complexity of the 
configuration, the results willhelp determine the validity of previous studies conducted at roomtemperature. 


II. Experimental Setup 
The experimental setup has four main elements: the rotating arm unit, the airfoil attached at the end of the am, 
the monosize droplet generator and the highs peed imaging s ystem. In Figure 1, the conceptual view of the experiment 
setup with all the elements except the high speed imaging system and the illumination is shown. The experimental 
setup components in the test cell before a run is presented in Figure 2. 


DROP GENERATOR 
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Figure 1: Conceptual view of the experiment. 
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Figure 2: Experimental setupin AERTS facility. 


A. The Adverse Environment Rotor Test Stand (AERTS ) 

The Adverse Environment Rotor Test Stand (AERTS) facility is located at Penn State University*. The facility 
was designed, andis used to test and evaluate rotor blade ice protections ystems. The AERTS facility is also used for 
ice protection coating evaluation and experimental rotor ice accretion shape correlation to ice shape modeling 
techniques. The facility is formed by an industrial6 m x 6 m x 6 m (20 ft x 20 ft x 20 ft) cold chamber where 
temperatures between -25°C and 0°C (-13°F to 32°F) can be achieved’. The chamber is cooled by convection using 
cooling lines and a set of fans located inside the chamber. Inside the chamber and surrounding the rotor there is a 
ballistic wall in the shape of an octagon. The ballistic wallis formed by 15.2 cm (6 in) thick weather resistant lumber 
reinforced with 0.635 cm (0.25 in) thick steel and covered by alummum plating for weather protection. A 93.2 kilo watt 
(125 hp) motor rotates the hub and the rotors. Figure 3 shows the chamber and its main components as seen from 
above. 


Figure 3: Adverse Environment Rotor Test Stand. 
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For the droplet breakup experiment, the AERTS facility was modified to accommodate the high speed vis ualization 
of droplet deformation and impact with the airfoil. The imaging systemwas located behind one of the chamber walls. 
An opening in the chamber wall allowed the cameras to observe and capture the droplet deformation and impact with 
the airfoil. A white 9500 lumens High Power LED from Cree was used forillumination. The LED was mounted on 
the hub of the rotor hub (Figure 4). The airfoils were located at the tip of the rotor blades. The rotor blades were 
truncated QH-S0 fiber-glass blades. The radius of the rotor was 1.57 m (62 in). During the experiment, the rotor 
operated from0 to 501 rpm. At the higher rotational s peed, the airfoil was traveling at 80 m/sec. 


Hub 


LED 
Mounted 
on Hub 


Figure 4: LED mounted on hub 


B. Airfoil Model 
The airfoil models mounted at the end of the rotor blades is a generic type of thick airfoil designated as DBKUP 
002 (Figure 5). Table 1 lists the coordinates of the airfoil. The airfoil was chosen forthe experiment because it has a 
blunt shape geometry that simulates a scaled version of the type of leading edge shape found on large transport airfoils. 
The airfoil measures 0.305 meters (12 inches ) in the spanwise direction, and has achord length of 0.47 meters (18.5 
inches ). 


Figure 5: DBKUP 002 airfoil model mounted on QH-50 rotor blade. 
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C. Monosize Droplet Generator 

A TSIMDG-100 Monosize Droplet Generator is used to produce water droplets within a diameter ran ge from 100 
um to 500 um. The droplet generator (Figure 6) is part of a system containing a water pressure container, a flow 
control valve, a frequency generator and the vibration head. Air from a compressor pressurizes the container and 
generates the water flow fromthe tank to the vibration head. A high precision flow rate control valve with a manoneter 
allows fine adjustment of the flow rate before the water flow reaches the vibration head. At the vibration head the 
water is forced through a small orifice (200 um in diameter) and a jet is formed. A disturbance in the formof a square 
wave at the appropriate frequency is introduced by activating a piezoelectric transducer controlled with a BK Precision 
model 4011A frequency generator. The jet is unstable at res onant frequencies and breaks into untformdroplets. Fora 
given orifice diameter, flow rate, and excitation frequency the diameter of the droplets generated is given by equation 
I: 


> 1/3 
D(um) = 317 (1) 
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Generator 
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Generator 
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Figure 6: TSI MDG-100 Monosize Droplet Generator. 


Where Qis the water flow rate in cubic centimeters per mmute; fis the frequency in kilohertz; and D1s the diameter 
of the droplets in micrometers. During the experiment orifices of 200 and 400um were used. The pressurized container 
consists of a5 gallon bucket with an airtight lid. The bucket was insulated and modified to accommodate a bendable 
immersion heater as wellas the 3.18 mm (1/8" in) OD air hose and tubing for the 12.7 mm (% in) OD, 6.35 mm (% in) 
ID water hose. Theimmersion heater was controlled by a PID controller to power on and off as necessary to maintain 
the desired water reservoir temperature. 


D. High Speed Imaging System 

The high speed imaging systemconsists of two high-speed cameras, the camera software, a triangular mirror and 
the lens system (Figure 7). The high-speed cameras used during the experiment were Photron SA-Z. The camera can 
capture images at rates from 1,000 to 1,000,000 frames per second (fps). During the experiment a frames rates of 
90,000 frames per second was used. A resolution of 384Hx492V pixels was used. The dimensions of the field of view 
were 11.45 millimeters in the horizontal directionby 14.66 millimeters in the vertical direction. It corresponded to a 
res olution of 29.8 micrometers per pixel. The camera s oftware employed for capturing the high speed image sequences 
is Photron’s FASTCAM Viewer (PFV) and is part of the camera system. The same software was used for post- 
processing the image sequences to generate the data movies for data analysis. During the experiment one lens 
configuration on each camera was used for magnification of the droplet deformation and/or breakup. The lens 
configuration consisted of a 200mm Micro Nikkor lens with a 2x doubler added between the camera andthe 200 mm 
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lens to double thefocallength ofthe microlens. Both cameras were placed 0.3m (11.88 inches) fromthe mirror. The 
distance fromthe tip of the mirror to the focal plane on the airfoil was 0.7 m (27.63 inches) (Figure 8). The lightning 
was designed to illuminate the droplets from behind to create a shadowgraph (black color of the droplets against a 
white-gray background). The light source was a white 9500 lumens High Power LED from Cree mounted on aring 


that was bolted to the blade grips. 


High Speed Cameras 
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Figure 7: Camera,lens, and triangular mirror setup. 
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Figure 8: Schematic of top view of cameraand triangular mirror placementrelative to airfoil model and droplet. 
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Ill. Test Procedure and Test Matrix 


A. Test Procedure 

Before a set of testruns, the high-speed camera s ystem was aligned, so that the line of view of the lens of each 
camera 1s parallel with the direction along the span of the airfoil at the midpoint of the leading edge (stagnation line 
along the span). Therotating armis set at the position, where the falling droplets fromthe monosize droplet generator 
will just graze the leading edge. The airfoilis held at the same location during the alignment process. The cameras are 
placed, with res pect to the trian gular mirror, at the distance required by the optics to obtain the needed magnification 
and such thatthe center of the lens depth of field is at the location where the falling droplets graze the airfoil leading 
edge. The motion of the cameras can be carefully controlled because each camera sits ona steel beamand on an x-y 
positioning table. 


The airfoil is moved from the alignment position and a ruler with the smallest subdivision of one millimeter 1s 
placed at the location where the droplets will graze the leading edgeof the airfoil. The camera is focused on theruler 
and the number of pixels per millimeter is recorded, sincethe number of pixels in the camera sensor and the field of 
view are known. This conversion is used during the data analysis to measure the droplet diameter. 


The monosize droplet generator controls are adjusted to produce the range of droplet sizes needed. The pressure 
and the frequency values are read fromthe calibration table obtained during the previous calibration of the s ystem for 
the orifice being used in the generator. The pressure valve and the frequency generator are set at the corresponding 
values andthe generatoris started. 


The pers onnel withdraw to the control room, where they are able to observethe motion and operation of the rotor 
through cameras. The setup cameras are controlled with software run on a laptop computer located in the control 
room. The software is used to start the camera and set the controlling parameters including the frame rate and the 
shutter speed. The software allows writing information that will be part of the camera frames. It is located on each 
frame above the recorded image. This information is helpful during the data analysis. During the experiment the 
following information was recorded on each frame: frame rate, time of recording, image resolution, frame number, 
lens configuration, date and time, airfoil velocity and target droplet size. 


The airfoils are set in motion at the revolutions per minute (rpm) corres ponding to the target velocity. When tt 
reaches the target speed, the camera recording is started and maintained for several rotations. Once the recording 1s 
completed therotoris brought to a stop. The high-speed movie is analyzed with the camera s oftware to determine the 
quality of the recording in each of the passes of the airfoil. If the recording is deemed of good quality the experiment 
moves to the next test point andthe process is repeated. 


B. Test Matrix 

The test matrix for the experiment is shown in Table 2. The test matrix was designed to obtain data for rotor tip 
speeds of 50, 60, 70, and 80 m/s and for droplet generator orifice sizes of 200 and 400 um. The pressure in the reservor 
and the frequency of the piezoelectric transducer were varied to obtain the desired droplet sizes when using the 
different orifices. The higher frame rates of 150000 and 140000 fps were used initially and compared to the lower 
frame rate of 90000 fps to determine the widest view available with a resolution that allowed for observation of the 
droplet from no deformation to deformation and breakup. The chamber temperature was maintained above freezng 
while adjusting the cameras to es tablish the testing setup. Once the cameras were ready to begin testing, the chamber 
temperature was decreased to create the environmentthat would allow the droplets to supercool at the impact location. 


IV. Droplet Supercooling Investigation 
Prior to collecting the high-speed video in each test, the temperature of the water is recorded as it exits the reservor 
and as it enters the droplet generator. A Type T thermocouple is inserted into the flow at the exit of the reservor and 
used with a PID controller connected to the immersed water heater to maintain the water at a constant desired 
temperature. The water in the reservoir needs to be maintained around 9°C to avoid freezing in the water hose and the 
droplet generator. A thermistor is inserted into the flow at the input of the droplet generator to measure the water 
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temperature as close as possible to the droplet jetting location. During the tests, the temperature of the water, after the 
droplets exited the generator, was not recorded to avoid interfering with the airfoil models. 


After the exact location of the droplet generator used in the experiments was established, the setup was repeated 
without spinning therotor to measure the temperature of the water at four vertical distances fromthe tip of the droplet 
generator. The temperature of the water was recorded at 2.54, 7.62, 12.7, and 19.05 cm (1, 3, 5, and 7.5 inches ) below 
the droplet generator using an infrared (IR) camera and a black surface. The impact of the water droplets against the 
airfoil occurred at 19.05 cm (7.5 in) below the jetting location. The resolution of the IR camera could not capture the 
droplets or streamto give a temperature reading. Therefore, a black surface, which was at the chamber temperature of 
-20°C, was used to allow the stream to impact on it allowing the IR camera to provide a temperature reading of the 
impacting water. The measurement approach estimates the temperature of the water droplets, but since heat transfer 
and convective cooling is expected to occur, uncertainty in the measurement exists. In the time between the fist 
contactof the water with the surface and the capturing of the temperature with the camera transient temperature effects 
occur, since the measurement does not measure the droplet temperature, butthat of the black surface subjected to the 
warming effect of the supercooled droplets. The mentioned shortcomings related to the measurement of the water 
droplet temperature using an infrared camera limits the validity of the datarecorded. Suchshortcomings are apparent 
by the large standard deviation observed on the experimental data, and the unexplained increase in droplet temperature 
at the 19.05 cm (7.5 in) location. For this reason, a prediction code® by Tadas P. Bartkus, a NASA Glenn Research 
Center contractor, was used to further confirm that the droplets were indeed supercooled. The relationship between 
the droplet temperature in degrees Celsius and the distance from the tip of the droplet generator as indicated by the 
experimental data and the predictions is dis played in Figure 9. The experimental data suggests that the temperature of 
the water is below freezing at each point below the tip. The predictions indicate that, while the droplet temperature 
may not be below freezing | inch below the tip, the water droplet temperature decreases to below O°C at the point of 
impact at 19.05 cm (7.5 inches ) below the tip of the droplet generator for droplets up to almost 700um in diameter. 
The diameter of all of the droplets observed was below 700um and therefore the droplets used in the experiment ate 
supercooled. 


The temperature at the tip of the droplet generator of 1°C was estimated as the initial droplet temperature in the 
models because the temperature using the thermistor, located 14.6 cm (5.75 in) above the tip, read between 2.65°C 
and 4.12°C. When the temperature at the entry of the droplet generator dipped below 2.65°C, the droplet generator tip 
would freeze. The temperature of the water 1s expected to decrease as it travels through the tubing adaptors and the 
droplet generator which are exposed to the chamber temperature, and given the freezing of the tip observed, the 1°C 
jetting temperature is considered representative. 


American Instituteof Aeronautics and Astronautics 


Experiment Twater=4C 
Simulation: 300mcr 


Experiment Twater=3C 
Simulation: 200mcr 


A A 

Bi Oo 

@ Simulation: 400mcr © Simulation: 500mcr 
@ Simulation: 600mcr O 


Simulation: 700 mcr 


Drop 73 


(°C) 


Distance from Droplet Generator 


(in) 


Figure 9: Droplet temperature at various distances from tip of droplet generator. For simulation purposes, droplet 
temperature at tip was approximated to be 1°C, with droplets maintaining a constant vertical velocity of 4 m/s. Droplet 
diameters range from 200 to 700 um. Chamber temperature is -20°C. 


V. Data Analysis 


A. Tracking aSingle Droplet to obtain the Droplet Dis placement against Time 

The data analysis is performed using the software program Photron FASTCAM Viewer 3 that allows frame by 
frame study of the movies, zooming in and out on single frames, running the movies forward and backward, and 
tracking a single droplet in x-y directions. Using PFV, a high speed movie corresponding to a set of experimental 
conditions is opened andthe frames where the airfoil impacts the droplets are located. There are ten sets of frames per 
movie that show the airfoil passage, because during thetest, the camera 1s programmed to capture ten rotations of the 
airfoil. In general several droplets can be observed impacting the airfoil and, in some cases, breaking up before impact, 
as sShownin Figure 10. 


Airfoil 


—| Tracked 
Droplet 


Figure 10: Droplet breakup before impact. Chamber temperature = -20°C, rotor tip velocity = 80 m/s, droplet diameter 
(left) =473 um, diameter (right) = 401 um. 
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The movie is run backward and forward to study the behavior of the droplets before and during their interaction 
with the airfoil. A droplet is selected for study and the programis returned to the frame where the airfoil is about to 
hit the droplet. The droplet may be deforming oris in the process of breaking up. A copy of the frame 1s saved in the 
experimental notes andthedroplet to be studied is labeled for future identification. The movie is run backward to the 
frame where there is no indication of deformation or interaction of the droplet with the approaching airfoil and where 
the tracking is going to be initiated, an example of both frames is shown in Figure 11. The frames that contain the 
motion of the droplet and the interaction with the airfoil are separated fromthe larger set of frames to create a shorter 
video. This shorter video is saved in .avifile format to be used and processed by aMATLAB program. 


) 
Airfoil moving 
from left to Tracked 
right. The ) Droplet 
droplet hasa 
vertical velocity 
of 3.84 m/s. 


Airfoil 


| 
Figure 11: First and last frame of video of a droplet being tracked. (Left: droplet tracking begins, Right: droplet tracking 
ends) Chamber temperature = -20°C, rotor tip velocity = 80 m/s, droplet diameter = 473 um, 165 framesin video, 1.83 
milliseconds in length. 


Once the file is in the correct format andin the file with the MATLAB script, the programrequires the frame rate, 
the number of fames in the video, the droplet number in the first frame (numbering fromtop to bottomthen going left 
to right, in the case of Figure 11, the droplet to be tracked is droplet number 3 in the frame), the number of the frame 
where the airfoil first appears, the airfoil velocity, and the resolution to convert from pixels to micrometers. The 
program begins by verifying the frame numbers and then calculating the appropriate threshold value to convert the 
grayscale images to binary images using the Otsu method’ and prepares the images for the subsequent parts in 
processing. The programthen determines the number of objects in the images and places an ellipse around the droplet 
being tracked. Thecentroid of the ellipse is also calculated and used to determine the motion of the droplet. From the 
tracking of the centroid, the motion in pixels relative to the first frame 1s recorded. For each frame, the frame number, 
the time with res pect to the first frame and the tracking frame, the horizontal and vertical dis placements of the centroid, 
the area, diameter, and perimeter of the ellipse, and the ellipse major and minor axes lengths are recorded. The 
information is then usedto calculate the droplet velocity and acceleration and the parameters of interest described in 
the following sections. 


B. Calculation of the Horizontal Velocity and Acceleration of the Droplet against Time 
A curve fit of the droplet dis placement against time is generated and used to calculate the horizontal velocity and 
acceleration of the droplet. The datapoints in the plot are curve-fit with a double five parameter exponential growth 
equation. The first and second derivatives of the curve-fit equation give the droplet velocity and acceleration. Those 
are the droplet horizontal velocity and acceleration against time measured by an observer located on a frame of 
reference at rest with res pect to the test cell. 


C. Change from the Frame of Reference at the Test Cell to the Frame of Reference at the Airfoil 
The horizontal dis placement, velocity, and acceleration were initially measuredin a frame of reference at rest with 
respectto thetestcell. It is more convenient to work in a frame of reference atrest on the airfoil. The frame of reference 
located on the airfoil has the origin located at the stagnation point on the leading edge. In a strict sense, a frame of 
coordinates at rest with respect to the airfoil is not an inertial frame because the airfoil is rotating, it is accelerated. 
Since the interaction between the droplets and the airfoil occurs at a distance less than or equal to the airfoil chord, 
that section of the airfoil path is nearly straight and the frame of reference at rest with respectto the airfoilis assumed 
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as inertial during that part of its path. The Galilean transformation together with the airfoil velocity was used to convert 


the velocity andacceleration of the droplet fromthe frame of reference at the test cell to the frame of reference at the 
airfoil (Figure 12). 
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Figure 12: Droplethorizontal velocity (top) and acceleration (bottom) in the frame of reference at reston the airfoil. 
Chamber temperature = -20°C, initial droplet velocity = 80 m/s, droplet diameter = 473 um. 


D. Position of a Droplet with respect to the Frame of Reference at the Airfoil 

Ata given location along the droplet path, the number of camera frames separating the droplet fromthe airfoil 1s 
equal to the frame number at which the airfoil hits the droplet mmus the frame number at the given location of the 
droplet. Multiplying the number of frames separating the droplet and the airfoil by the time lapse between frames 
(1/frame rate) and by the airfoils peed gives the horizontal position of the droplet with respect to the frame of reference 
at the airfoil. Repeating thecalculation for each droplet location gives the distance of the droplet along its horizontal 
path with respect to the leading edge of the airfoil. The main parameters in the experiment were calculated at each 
position of the droplet as it approaches the airfoil and with respect to the frame of reference located at rest on the 
airfoil. 
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E. Air Velocity at a given Droplet Location 

For an airfoil with the same geometry as the DBKUP 02 but scaled down 1/10 in size (0.047m chord), the air 
velocity along the horizontal line ending at the stagnation point on the leading edge was measured experimentally 
with a Laser Droplet Velocimeter (LDV) at the INTA 0.3mx0.2m wind tunnel. The vertical axis in the plot records 
the air velocity divided by the free streamvelocity (Vair/ Va). The horizontal axis records the distance with respect to 
the leading edge of the airfoil divided by the chord (x/c). The experimental data was curve-fit and the resulting equation 
was used to calculate the air velocity at any given droplet position along its path. To calculate the position of the 
droplet using the curve-fit equation the value of the free stream velocity is needed. For the frame of reference at test 
with respectto the airfoil, the air free streamvelocity is the same as the velocity of the airfoil with res pect to the frame 
of reference at rest with respectto the test cell. 


F. Relative Velocity between the Droplet and the Air (Slip Velocity) 

Since at each droplet position along its path the horizontal velocity of the droplet and the velocity of the air are 
known, the relative velocity between the droplet and the air can be calculated. This velocity, | Vip Va noniel | is called 
the "slip velocity”. The three velocities, the droplet, air, and slip velocities, are shown in Figure 13. The Reynolds 
number and the Weber number can be calculated when the horizontal slip velocity at each location along the dropket 
path is known. 
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Figure 13: Droplet, air, and slip Velocity against distance from the leading edge of the airfoil. Chamber temperature = 
-20°C, initial droplet velocity =80 m/s, droplet diameter = 473 um. 


G. Calculation of the Reynolds Number, the Weber Number 

The main parameters measured along the path of the droplet were the Reynolds number (Eqn. 2), the Weber 
number (Eqn. 3), and the Bond Number (Eqn. 4). All those parameters were measured with res pect to the horizontal 
motion of the droplet. The droplets leave the monosize droplet generator and fall vertically at a constant speed between 
3-4 m/sec (measured during the experiment). There was no indication fromthe physical behavior of the droplet that 
the vertical velocity had an influence on the deformation of the droplet. The shape of the droplet during deformation 
corres ponded to the response of the droplet to the horizontal velocity and acceleration. In all the measurements and 
calculations the vertical velocity was not considered. 


The following are the definitions of the main parameters measured along the path of the droplet: 


R = Pair Van Varina D 
Reynolds Number or Hair (2) 
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Weber Number We OREN ee (3) 
2( dV 
pie Pied) droplet 
Bond Numb O water / air d 
ond Number f (4) 


Where Dis the droplet diameter, Varopieris the droplet velocity, Vair is the air velocity at the location of the dropket, 
(Vair — Varopiet 18 the slip velocity, pair is the air density, fair is the air absolute viscosity, Pwarer 1S the water density for 
the droplet, Owaternir 18 the water surface tension for the droplet. The Bond number was defined with respect to the 
droplet acceleration. 


Since the air and water droplet properties are known and the droplet diameter, the slip velocity and the droplet 
acceleration are measured or calculated along the path of the droplet, then the non-dimensional parameters can be 
calculated. The velocities and accelerations employed in the calculation of the parameters are for the horizontal 
direction. 


VI. Results: Comparison of Supercooled and Room Temperature Droplets 


A. Comparison of Supercooled Droplet and Room Temperature Droplet Data Collectin the AERTS Facility 
and in INTA 

For the comparison of the droplet behavior, data fromthe tests done at INTA at roomtemperature are compared 
to data from the tests done in the AERTS facility at -20°C chamber temperature. Although the airfoil speeds for the 
tests are different, 90 m/s for INTA tests and 80m/s forthe AERTS facility tests, the behavior of the droplets ata s lip 
velocity of 60 m/s is compared. The deformation (ellipse major axis/ellipse mmor axis) and the Bond number of the 
droplets are plotted against the droplet diameter in Figure 14. The tests at INTA used droplet generator orifice sizes 
ranging from 500 to 1000 um while 200, and 400 um droplet generator orifice sizes were used when testing in the 
AERTS facility. The lack of overlap in the datadoes not allow for proper comparis ons to be made between the dropkts 
in relation to their behavior. Yet, both the supercooled droplets and the roomtemperature droplets behave in a similar 
manner in that larger diameter droplets have higher deformation than smaller droplets. This can be observed via visual 
inspection. Two distinct sets of data related to the warm temperature case and the supercooled case arise in the 
comparison presented in Figure 14, delimited by the line sketched in the graph. 


Given that the described comparisons are conducted at varying airfoil velocities, environmental temperatures, and 
droplet sizes, a conclusion regarding the effects of droplet supercooling on the deformation of the drops cannot be 
concluded. However, let it be assumed that the supercooling of the droplets has no effect on the droplet fragmentation. 
If this is thecase, the Bondnumber would not change. Water surface tension is the main parameter in the Bond number 
that would vary with temperature. The assumption of equal droplet deformation behavior for both environmental 
conditions assumes thatthe surface tension of the drops can be estimated to be constant giventhesmall (less than 5%) 
variation of water surface tension between room temperature and -10°C water droplets. If this is the case, the Bond 
number can be plotted versus droplet size (Figure 15). The results suggests that there 1s a velocity effect between the 
two sets of data investigated (room temperature vs. supercooled droplets). To further explore the reason for the 
difference between the two environmental conditions explored, the deformation of droplets at the same temperature 
for varying airfoil velocities (initial droplet velocity in the fixed frame) are compared in the next section. 
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trendlines were added to assist with the interpretation of the results. The slopes obtained are in agreement in 
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Figure 14: Droplet deformation against diameter at a slip velocity of 60 m/s. 
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Figure 15: Bond number against droplet diameter at a slip velocity of 60 m/s. 


B. Effect of Initial Velocity on Deformation and Bond Number 

The deformation for the impacting drops tested at varying airfoil velocities plotted against the droplet diameter are 
shown in Figure 16. For these cases, the air temperature is -20°C and the slip velocity is 50 m/s, but the initial droplet 
velocities related to the airfoil velocity in the rotating frame are 80 m/s, 70 m/s, and 60 m/s. Such comparison allows 
for the investigation of theeffects of velocity (initial and slip) on the droplet deformation. In the conducted correlation, 


magnitude, and there is no indication of a particular trend with variations in velocity (See Figure 16). 


To investigate the effects of slip velocity, the deformation against droplet diameter for droplets with initial velocity 
of 70 m/s and at fourslip velocities at a chamber temperature of -20°C is plotted in Figure 17. Foran initial droplet 


velocity of 70 m/s, the droplet deformation slope at a slip velocity of 30m/s is 0.0011. As the slip velocity of the drop 
increases to 60 m/s, the droplet deformation slope increases to 0.0065. The same effect on the droplet deformation 


slope is observed as the slip velocity approaches the airfoil velocity. These results indicate that the ratio of the slip 
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velocity to the initial velocity affects the slope of the droplet deformation with res pect to the droplet size. Given such 
variations in the droplet deformation slopes with respect to the droplet size, the initial droplet velocity and the skp 
velocity must be the same to investigate the effects of supercooling on droplet deformation. The following section 
compares the deformation of droplets with an initial drop velocity of 60.m/s at a slip velocity of 40 m/s at two varying 
environmental temperatures in an attempt to determine theeffects of water temperature on such droplet deformation. 
The selection of mentioned slip and initial velocities ensures an equal ratio of these velocities for the investigation of 
the effects of supercooling on droplet deformation. 
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Figure 16: Dropletdeformation against diameter at a slip velocity of 50 m/s and chamber temperature of -20°C. 
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Figure 17: Deformation against droplet diameter for droplets with initial velocity of 70 m/s and at four slip velocities, 
chamber temperature = -20°C. 
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C. Comparison of Supercooled Droplets and Room Temperature Droplet Data Collected in the 
AERTS Facility 

It has been established thatto determine the effects of droplets upercooling onits deformation prior to impact with 
an airfoil, identical slip and droplet initial velocities mustbe compared. In addition, an effort is made to also compare 
data containing similar droplet diameters. The datacollected using an airfoil speed of 60 m/s in the AERTS facility at 
20 “C (room temperature) and -20°C were used. In the mentioned tests the orifices used in the droplet generator were 
150 and 200 micrometer diameters. The droplet deformation at a slip velocity of 40 m/s is plotted below in Figure 18. 
The data in the figures overlaps as all of the droplets analyzed have similar diameters and trends indicating differences 
between warmdroplets and supercooled droplets is not observed for the tested conditions. Based on theseresults, it is 
expected that the Bond number for the two environmental conditions should also agree, suggesting that the variation 
of the surface tension of the droplets does not cause a change in the behavior of the supercooled droplets. For this 
reason, it is assumed that the surface tension of the water is constant for both environmental conditions. If this is the 
case, itis expected thatthe Bondnumber for both temperature cases should follow an equal trend. The bond number 
for the room temperature case and the supercooled environment is plotted in Figure 19. Both the warm environment 
and the supercooled case follow an equal trend. Square polynomnal fittings (since Bond number is proportional to the 
square of the droplet diameter) were added, further showcasing the continuous trend between the results. The 
processed experimental datas uggests that supercooling has no effect on the particle deformation for the airfoil velocity 
mves tigated. 
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Figure 18: Dropletdeformation against diameter at a slip velocity of 40 m/s and initial droplet velocity of 60 m/s. 
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Figure 19: Bond number against diameter at a slip velocity of 40 m/s and initial droplet velocity of 60 m/s. 


VII. Conclusions 

The purpose of the work presented here was to determine the effects of supercooling on the behavior of water 
droplets near the leading edge of an airfoil. The room temperature droplet data that was collected at INTA in Spain 
was compared to the data collected of supercooled droplets in the AERTS facility. The two sets of data could not be 
used in the comparison of the behavior as the data did not overlap in terms of slip velocity/initial velocity ratios. 
Although thesame slip velocity was compared, theranges of droplet diameters and the airfoil velocity were different. 
It was demonstrated that to compare the effects of water supercooling on droplet deformation, the ratio of the slip 
velocity and the initial droplet velocity must be equal. To further explore the effects of droplet velocity on the 
deformation of droplets, experimental cases with constant slip velocity and varying airfoil velocity were investigated. 
A specific trend was notobservedin the slopes of the data, but when the airfoil velocity was held constant and the shp 
velocity was varied, a clear effect on droplet deformation was observed. As the slip velocity approached the airfoil 
velocity, the slope of the data increased. For an airfoil velocity of 70 m/s and a slip velocity of 30 m/s the slope was 
0.0011 and it increased to 0.0065 at a slip velocity of 60 m/s. This indicated that the ratio of the slip velocity to the 
airfoil velocity should be the same when comparing datato determine the effects of supercooling on droplet behavior. 


A case with equal slip velocity to initial velocity ratios was selected for roomtemperature and supercooled dropkt 
conditions. The airfoil velocity was 60m/s and the slip velocity for both sets of data was 40 m/s. In these cases, the 
deformation of the supercooled droplet and that of the roomtemperature droplets did not present differenttrends for 
the limited range, weakly supercooled states of water droplets, tested. The similar behavior for both environmental 
conditions indicates that weak supercooling had no effect on particle deformation at the selected impact velocity. The 
assumption of a constant surface tension value for the mentioned conditions was further supported by the equal trend 
of the Bond number obtained for supercooled droplets and roomtemperature ones. 


Acknowledgments 


The authors would like to thank Mr. Chris Lynch and Mr. Quentin Schwinn for their excellent work with the high- 
speed imaging system. We would also like to thankall of the AERTS facility research assistants that were involved 
in the setup of the experiments and collection of data. Your help is very much appreciated. 


18 
American Instituteof Aeronautics and Astronautics 


References 


' Tan, J., Papadakis, M.,and Sampaht, M. K., “Computational Study of Large Droplet Breakup in the Vicinity of 
an Airfoil”, DOT/FA A/AR-05/42, Final Report, October 2005. 

*Feo, Alex., Vargas, M., and Sor, S., “Rotating Rig Development for Droplet Deformation/Breakup and Impact 
Experiments Induced by Aerodynamic Surfaces”, to be published in 2011. 

*Vargas, M., Sor, S., and Magarino, A.G., “Mechanismof Water Droplet Breakup near the Leading Edge ofan 
Atrfoil,’ AIAA Paper 2012-3129, June 2012. 

“Soltis, J., Palacios, J., Eden, T., and Wolfe, D., “Ice Adhesion Mechanisms of Erosion -Resis tant Coatings,” AIAA 
Journal, Early Edition, June 19, 2014. 

doi: 10.2514/1.J053208 

*Palacios,J., Han, Y., Brouwers, E. W., Smith, E. C., “Icing EnvironmentRotor Test Stand Liquid W ater Content 
Measurement Procedures and Ice Shape Correlation,” Journal ofthe American Helicopter Society, Vol. 57, (2), April 
2012, pp. 29-40. 

° Bartkus, T. P., Struk, P. M., Tsao, J. and Van Zante, J. F., “Numerical Analysis of Mixed -Phase Icing Cloud 
Simulations in the NASA Propulsion Systems Laboratory,” 8th AIAA Atmospheric and Space Environments 
Conference, AIAA, Washington D.C., 2016, AIAA 2016-3739. 

7Otsu,N.,"A Threshold Selection Method from Gray -Level Histograms," IEEE Transactions on Systems, Man, 
and Cybernetics, Vol.9, No. 1, 1979,pp. 62-66. 


19 
American Instituteof Aeronautics and Astronautics 


X-Coordinate 


1.0000 
0.9989 
0.9957 
0.9904 
0.9830 
0.9735 
0.9619 
0.9484 
0.9330 
0.9157 
0.8967 
0.8759 
0.8536 
0.8297 
0.8044 
0.7778 
0.7500 
0.7211 
0.6913 
0.6607 
0.6294 
0.5976 
0.5653 
0.5327 
0.5000 
0.4673 
0.4347 
0.4025 
0.3706 
0.3393 
0.3087 
0.2789 
0.2500 
0.2222 
0.1956 
0.1703 
0.1465 
0.1241 
0.1033 
0.0843 
0.0670 
0.0516 
0.0381 
0.0265 
0.0170 
0.0096 
0.0043 
0.0011 


Y-Coordinate 


0.0000 
-0.0009 
-0.0014 
-0.0025 
-0.0042 
-0.0064 
-0.0089 
-0.0120 
-0.0155 
-0.0195 
-0.0240 
-0.0291 
-0.0347 
-0.0410 
-0.0479 
-0.0554 
-0.0634 
-0.0721 
-0.0815 
-0.0915 
-0.1019 
-0.1128 
-0.1240 
-0.1353 
-0.1466 
-0.1577 
-0.1683 
-0.1780 
-0.1870 
-0.1938 
-0.1985 
-0.2006 
-0.2007 
-0.1985 
-0.1943 
-0.1880 
-0.1797 
-0.1700 
-0.1584 
-0.1458 
-0.1318 
-0.1170 
-0.1007 
-0.0842 
-0.0674 
-0.0498 
-0.0310 
-0.0139 
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Table 1. DBKUP 002 Airfoil Coordinates 


X-Coordinate 


0.0000 
0.0011 
0.0043 
0.0096 
0.0170 
0.0265 
0.0381 
0.0516 
0.0670 
0.0843 
0.1033 
0.1241 
0.1465 
0.1703 
0.1956 
0.2222 
0.2500 
0.2789 
0.3087 
0.3393 
0.3706 
0.4025 
0.4347 
0.4673 
0.5000 
0.5327 
0.5653 
0.5976 
0.6294 
0.6607 
0.6913 
0.7211 
0.7500 
0.7778 
0.8044 
0.8297 
0.8536 
0.8759 
0.8967 
0.9157 
0.9330 
0.9484 
0.9619 
0.9735 
0.9830 
0.9904 
0.9957 
0.9989 


Y-Coordinate 


0.0000 
0.0139 
0.0310 
0.0498 
0.0674 
0.0842 
0.1007 
0.1170 
0.1318 
0.1458 
0.1584 
0.17700 
0.1797 
0.1880 
0.1943 
0.1985 
0.2007 
0.2006 
0.1985 
0.1938 
0.1870 
0.1780 
0.1683 
0.1577 
0.1466 
0.1353 
0.1240 
0.1128 
0.1019 
0.0915 
0.0815 
0.0721 
0.0634 
0.0554 
0.0479 
0.0410 
0.0347 
0.0291 
0.0240 
0.0195 
0.0155 
0.0120 
0.0089 
0.0064 
0.0042 
0.0025 
0.0014 
0.0009 


Table 2: Test Matrix 


ee Chamber 
Rotor Velocity Orifice |Pressure| Frequency | Thermistor 
Date Temperature | Temperature FPS 
os (m/s) Size (um) (°C) (°C) 


2/21/17 14.00 15.90 20.00 140000 
14.00 17.90 20.00 140000 
14.00 17.20 20.00 140000 
50.00 17.30 20.00 90000 
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